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The extended kernel ridge regression (EKRR) method with odd-even effects was adopted to improve the de- 
scription of the nuclear charge radius using five commonly used nuclear models. These are: (i) the isospin depen- 
dent A!/? formula (Gii) relativistic continuum Hartree-Bogoliubov (RCHB) theory (iii) Hartree-Fock-Bogoliubov 
(HFB) model HFB25 (iv) the Weizsicker-Skyrme (WS) model WS“, (v) HFB25* model. In the last two models, 
the charge radii were calculated using a five-parameter formula with the nuclear shell corrections and deforma- 
tions obtained from the WS and HFB25 models, respectively. For each model, the resultant root-mean-square 
deviation for the 1014 nuclei with proton number Z > 8 can be significantly reduced to 0.009-0.013 fm af- 
ter considering the modification with the EKRR method. The best among them was the RCHB model, with a 
root-mean-square deviation of 0.0092 fm. The extrapolation abilities of the KRR and EKRR methods for the 
neutron-rich region were examined and it was found that after considering the odd-even effects, the extrapo- 
lation power was improved compared with that of the original KRR method. The strong odd-even staggering 
of nuclear charge radii of Ca and Cu isotopes and the abrupt kinks across the neutron N = 126 and 82 shell 


closures were also calculated and could be reproduced quite well by calculations using the EKRR method. 
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I. INTRODUCTION 


The nuclear charge radius, similar to other quantities such 
as the binding energy and half-life, is one of the most basic 
properties reflecting the important characteristics of atomic 
nuclei. Assuming a constant saturation density inside the nu- 
cleus, the nuclear charge radius is usually described by the 
A‘/® law, where A is the mass number. By studying the 
charge radius, information on the nuclear shells and subshell 
structures [1, 2], shape transitions [3, 4], the neutron skin and 
halos [5-7], etc., can be obtained. 

With improvements in the experimental techniques and 
measurement methods, various approaches have been adopt- 
ed for measuring the nuclear charge radii [8, 9]. To date, 
more than 1000 nuclear charge radii have been measured 
[10, 11]. Recently, the charge radii of several very exot- 
ic nuclei have attracted interest, especially the strong odd- 
even staggering (OES) in some isotope chains and the abrupt 
kinks across neutron shell closures [2, 12-21], which provide 
a benchmark for nuclear models. 

Theoretically, except for phenomenological formulae [22— 
29], various methods, including local-relationship-based 
models [30-35], macroscopic-microscopic models [36-39], 
nonrelativistic [40—43] and relativistic mean-field mod- 
el [44-52] were used to systematically investigate nuclear 
charge radii. In addition, the ab-initio no-core shell mod- 
el was adopted for investigating this topic [53, 54]. Each 
model provides fairly good descriptions of the nuclear charge 
radii across the nuclear chart. However, with the exception 
of models based on local relationships, all of these methods 
have root-mean-square (RMS) deviations larger than 0.02 fm. 
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It should be noted that few of these models can reproduce 
strong OES and abrupt kinks across the neutron shell closure. 
To understand these nuclear phenomena, a more accurate de- 
scription of nuclear charge radii is required. 


Recently, due to the development of high- performance 
computing, machine learning methods have been widely 
adopted for investigating various aspects of nuclear physic- 
s [55-59]. Several machine learning methods have been used 
to improve the description of nuclear charge radii, such as 
artificial neural networks [60-63], Bayesian neural network- 
s [64-68], the radial basis function approach [69], the ker- 
nel ridge regression (KRR) [70], etc. By training a machine 
learning network using radius residuals, that is, the deviations 
between the experimental and calculated nuclear charge radii, 
machine learning methods can reduce the corresponding rms 
deviations to 0.01-0.02 fm. 


The KRR method is one of the most popular machine- 
learning approaches, with the extension of ridge regression 
for nonlinearity [71, 72]. It was improved by including odd- 
even effects and gradient kernel functions and provided suc- 
cessful descriptions of various aspects of nuclear physics, 
such as of the nuclear mass [73-77], nuclear energy den- 
sity functionals [78], and neutron-capture reaction cross- 
sections [79]. In the present study, the extended KRR (EKR- 
R) method with odd-even effects included through remodula- 
tion of the KRR kernel function [74] is used to improve the 
description of the nuclear charge radius. Compared with the 
KRR method, the number of weight parameters did not in- 
crease in the EKRR method. 


The remainder of this paper is organized as follows. A 
brief introduction to the EKRR method is presented in Sec. 
II. The numerical details of the study are presented in Sec. 
HI. The results obtained using the KRR and EKRR methods 
are presented in Sec. IV. The extrapolation power of the 
EKRR method is discussed. The strong OES of the nuclear 
charge radii in Ca and Cu isotopes and abrupt kinks across the 


67 neutrons N = 126 and 82 shell closures were investigated. 
e Finally, a summary is presented in Sec. V. 


69 Il. THEORETICAL FRAMEWORK 


77 The KRR method was successfully applied to improve 
7 the descriptions of nuclear charge radii obtained using 
72 several widely used phenomenological formulae [70]. 
73 To include odd-even effects, the KRR function S() 
u X. K(zj, xi)a;i is extended to the EKRR function [74] 


m m 
75 S(2x;) =) 0 K(2;,01)0; + > Koe(x;,21)Bi , Gd) 
i=1 i=1 
ze Where x; are the locations of the nuclei in the nuclear chart, 
7 with x; = (Zi, Ni). m is the number of training data points, 
7 a; and (; are the weights, K(x;,2;) and Koe(£j, £i) are k- 
79 ernel functions that characterize the similarity between the 
so data. In this study, a Gaussian kernel was adopted, which is 
a31 expressed as 


(2) 


ea where ||z; — x;|| = \/(Zi — Z;)? + (Ni — Nj)? is the dis- 
sa tance between two nuclei. Koe(X;, zi) was introduced to en- 
ss hance the correlations between nuclei with the same number 
ss parity of neutrons and protons, which can be written as: 


82 K (aj, x) = exp(—||ai — z;||*/20°) , 


(3) 


s Ôoe(£j, £i) = 1 (0) if the two nuclei have the same (differ- 
s ent) number parities of protons and neutrons. o and Coe are 
əx hyperparameters We defined the range affected by the kernel. 
ə% The kernel weights a; and 8; are determined by minimiz- 
ə2 ing the following loss function: 


a Koeltj, £i) = Soe(aj, 24)exp(—||24 — 25||?/203.) - 


L(a, B) = 5 [S(z:) — y(z:)]? + aT Ka + ocb T KoeB . 
i=1 

93 (4) 
ot The first term is the variance between the training data y(x;) 
ə and the EKRR prediction S(«;). The second and third terms 
æ are regularizers, where the hyperparameters \ and Ase deter- 
97 mine the regularization strength and are adopted to reduce the 
ə risk of overfitting. 
ə By minimizing the loss function [Eq. (4)], we obtain 


À 
100 B = A 5 (5) 
A —1 
101 a= (K+ Koo +41) yY. (6) 
102 According to Eq.(5), the EKRR function [Eq. (1)] can be 
103 Written as a standard KRR function: 
104 S(x;) = 5 K'(z;, zi)@i , (7) 


os where K'(x;, x;) is the remodulation kernel. 


(8) 
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107 According to Eq.(5), the number of weight parameters in the 
108 EKRR method is identical to that in the original KRR method. 


K' (xj, xi) = K (aj, xi) + 


09 II. NUMERICAL DETAILS 


uo In this study, 1014 experimental data points with Z > 8 
111 were considered and obtained from Refs. [10, 11]. The EKR- 
112 R function (7) was trained to reconstruct the residual radius: 
113 i.e., the deviations AR(N, Z) = R®®(N, Z) — R®(N, Z) 
114 between the experimental data R°P(N, Z) and the predic- 
us tions Rt? (N, Z) for the following five nuclear models. 


(i) The widely used phenomenological formula R, 
ra [1 — b(N — Z)/A] A!/3 [24] with the parameter 
ra=1.282 fm and b = 0.342 was fitted by experimen- 
tal data (further denoted by A!/*). 


(ii) The relativistic continuum Hartree-Bogoliubov (RCHB) 
theory [47]. 


(iii) The MHartree-Fock-Bogoliubov (HFB) model 
B25 [80]. 


HF- 


(iv) The Weizsacker-Skyrme (WS) model WS* [11]. 
(v) The HFB25* model [11]. 


Note that by considering the nuclear shell corrections and 
deformations obtained from the WS and HFB25 models, a 
five-parameter nuclear charge radii formula was proposed in 
Ref. [11]. In this study, these methods are denoted as WS* 
and HFB25*, respectively. The parameters in the formulae 
of these two models were obtained from Refs. [11]. The 
RMS deviations between the experimental data and the five 
models (A,ms) are listed in Table 1. Once the weights a; 
were obtained, the EKRR function S(N, Z) was obtained for 
each nucleus. Therefore, the predicted charge radius for a 
nucleus with neutron number NV and the proton number Z is 
given by REKRR = Rt(N, Z) + S(N,Z). In this study, 
the KRR method was adopted for predicting charge radii for 
comparison. 

Leave-one-out cross-validation was adopted to determine 
the two hyperparameters (o and A) in the KRR method and 
the four hyperparameters (o, A, Coe and Age) in the EKRR 
method. The predicted radius for each of the 1014 nuclei can 
144 be given by the KRR/EKRR method trained on all other 1013 
145 nuclei with a given set of hyperparameters. The optimized 
146 hyperparameters (see Table 1) are obtained when the RMS 
deviation between the experimental and calculated radii reach 
a minimum value. 


147 


IV. RESULTS AND DISCUSSION 


Table 1 lists the hyperparameters (a, A) in the KRR method 
and (a, A, Coe and Age) using the EKRR method as well as 
152 the RMS deviations between the experimental data and the 
153 predictions of the five models. The RMS deviations with 
154 (without) KRR and EKRR are denoted by AKER 
iss (Arms). With the exception of the phenomenological A!/3- 
156 formula all other models provided a good global description 
157 Of the nuclear charge radii, especially for the WS*. It should 
iss be noted that a spherical shape is considered in the RCHB 


TABLE 1. The hyperparameters (o, A, Coe and Ace) in the KRR and EKRR method, and the RMS deviations between the experimental 


data and the predictions by five different models. The RMS deviations with (without) KRR and EKRR methods are denoted by AFPR and 
A (Arms): 
Model o AÀ Joe doe Arms (fm) Ane (fm) A R (fm) 
Avs = à = £ 0.0672 - - 
2.84 0.01 - - - 0.0158 - 
2.32 0.01 2.88 0.02 - - 0.0100 
RCHB - - - - 0.0350 [47] - - 
2.68 0.02 - - - 0.0157 - 
1.83 0.01 2.73 0.02 - - 0.0092 
HFB25 - - - - 0.0256 [80] - - 
1.77 0.34 - - - 0.0177 - 
1.48 0.08 2.20 0.22 - - 0.0130 
WwSs* - - - - 0.0210 [11] - - 
0.70 0.01 - - - 0.0155 - 
1.54 0.02 2.46 0.03 - - 0.0096 
HFB25* - - - - 0.0254 [11] - - 
0.68 0.01 - - - 0.0182 - 
1.35 0.05 2.21 0.08 - - 0.0120 


159 theory when investigating the entire nuclear landscape [47]. 
160 Therefore, its RMS deviation is slightly larger than that for 
161 the nonrelativistic model HFB25. To date, only even-even nu- 
162 Clei have been calculated in the deformed relativistic Hartree- 
163 Bogoliubov theory on a continuum (DRHBc) [48, 51]. The 
164 description of the nuclear charge radii can be further im- 
165 proved when all nuclei in the nuclear chart are calculated 
166 using this model. It can also be observed that HFB25 and 
167 HFB25* yield similar RMS deviations when describing the 
nuclear charge radii. After the KRR method had been con- 
sidered, all RMS deviations for these five models could be 
significantly reduced to approximately 0.015-0.018 fm, par- 
ticularly for the A‘/* formula. Interestingly, the RMS devi- 
ations of the HFB25 and HFB25* models were smaller than 
73 those of the A!/* formula and the RCHB model without the 
17 KRR method. However, after the KRR method was consid- 
175 ered, the situation was reversed. After considering the odd- 
176 even effects, the predictive powers of the five models were 
177 further improved by the EKRR method compared with the 
172 KRR method. The RMS deviation was further reduced by ap- 
179 proximately 0.006 fm for the five models, with the exception 
180 Of the HFB25 model, for which it was reduced to less than 
181 0.005 fm. The RMS deviations of the three models (A!/* for- 
182 mula, RCHB and WS*) were less than 0.01 fm, whereby the 
183 smallest was for the RCHB model with an RMS deviation e- 
184 qual to 0.0092 fm. This is the best result for nuclear charge 
185 radii predictions using the machine learning approach, as far 
186 aS We are aware. Here, we show the typical RMS deviations 
of some popular machine learning approaches. 
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87 
188 (i) artificial neural network: 0.028 fm [61]; 


(ii) Bayesian neural network: 0.014 fm [68]; 


89 


(iii) radial basis function approach: 0.017 fm [69]. 


90 


ə1 Note that if the full nuclear landscape is calculated using the 
DRHBc theory, the description of the nuclear charge radii can 


still be improved using the EKRR method. To show Table 1 
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A? RCHB HFB25 WS* HFB25* 
Fig. 1. (Color online) The RMS deviations between the experimental 
data and the predictions of five different models with and without the 


KRR/EKRR method. 


194 in a more visual manner, a comparison of these five models is 
also shown in Fig. 1. 

Figure 2 shows the differences in the radii between the 
experimental data and the calculations of the RCHB mod- 
el (grey solid circles), KRR method (red triangles) and the 
19 EKRR methods (blue crosses). Because the improvements 
achieved by the KRR and EKRR methods for the five models 
mentioned above were similar, we consider only the RCH- 
B model as an example. In order to study the odd-even ef- 
fects included in the EKRR method, the data were divided 
into four groups characterized by even or odd proton num- 
205 bers Z and neutron numbers JN, that is, even-even, even-odd, 
odd-even, and odd-odd. Clearly, the predictive power of the 
RCHB model could be further improved by using the EKRR 
method compared with the original KRR method. The signif- 
icant improvement of the EKRR method is mainly due to the 
consideration of the odd-even effects, which eliminates the 
staggering behavior of radius deviations owing to the odd and 
even numbers of nucleons using the KRR method. It can be 
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Fig. 2. (Color online) Radius differences AR between the experi- 
mental data and the calculations of the RCHB model (grey solid cir- 
cles), the KRR method (red triangles), and the EKRR method (blue 
crosses) for (a) even-even, (b) even-odd, (c) odd-even, and (d) odd- 
odd nuclei. 


seen that when the mass number is A ~ 150, the prediction- 
s of the KRR method exhibit significant deviations from the 
data, which can be significantly improved using the EKRR 
method. This is clear evidence of the importance of consider- 
ing the odd-even effects in predictions of the nuclear charge 
218 radius. 

219 To investigate the extrapolation abilities of the KRR and 
22 EKRR methods for neutron-rich regions, the 1014 nuclei with 
221 known charge radii were redivided into one training set and 
six test sets as follows: For each isotopic chain with more 
than nine nuclei, the six most neutron-rich nuclei were se- 
lected and classified into six test sets based on the distance 
from the previous nucleus. Test set 1 (6) had the shortest 
(longest) extrapolation distance. This type of classification is 
the same as that used in our previous study [70]. The hy- 
223 perparameters obtained by leave-one-out cross-validation in 
the KRR/RKRR method remained the same in the following 
calculations: 

RMS deviations of the KRR and EKRR methods for dif- 
ferent extrapolation steps for the five models are shown in 
Figs. 3(a)-(e). A clearer comparison of the RMS deviations 
scaled to the corresponding RMS deviations of the five mod- 
els without KRR/EKRR corrections are shown in Figs. 3(f) 
and (j). Regardless of whether the KRR or EKRR method is 
considered, the RMS deviation increased with the extrapola- 
tion distance. For the A!/? formula and the RCHB model, 
the KRR/EKRR method could improve the radius descrip- 
tion for all extrapolation distances. For the other three mod- 
els, the KRR method only improved the radius description for 
an extrapolation distance of one or two, which could be fur- 
ther improved after considering the odd-even effects with the 
EKRR method. This indicates that the KRR/EKRR method 
loses its extrapolation power at extrapolation distances larg- 
er than 3 for these three models. This is due to the charge 
247 radii calculated using these three models, which were quite 
good, and their RMS deviations, which were already suffi- 
ciently small. The KRR/EKRR method automatically iden- 


248 


24! 


© 


250 tifies the extrapolation distance limit owing to the hyperpa- 
251 rameters o and Coe being optimized using the training data. 
22 Refs. [73, 74] demonstrated that the KRR and EKRR method- 
253 S lose their predictive power at larger extrapolation distances 
254 (approximately six), when predicting the nuclear mass using 
255 the mass model WS4 [81]. This may be due to much more 
256 mass data existed than the charge radii, and the KRR/EKRR 
257 networks can be trained better with more data. In general, the 
28 EKRR method has a better predictive power than the KRR 
method for an extrapolation distance of less than 3. For an 
extrapolation distance greater than 3, the results of the KRR 
and EKRR methods were similar in most cases. Almost none 
of these extrapolations exhibited overfitting, except for WS* 
at an extrapolation distance of 3, and this overfitting was quite 
small. This indicates that both the KRR and EKRR method- 
s have good extrapolation powers and can avoid the risk of 
overfitting to a large extent. 

The observation of the strong OES of the charge radii 
throughout the nuclear landscape provides a particularly strin- 
gent test for nuclear theory. To examine the predictive power 
of the EKRR method, which is improved by considering the 
odd-even effects compared with the original KRR method, in 
the following we will investigate the recently observed OES 
of the radii in calcium and copper isotopes [14—16]. Similar 
to the gap parameter, the OES parameter for the charge radii 
is defined as: 


259 


1 
A® (Z, N) = gi"(Z,N -1) - 2r(Z, N) +r(Z,N + 1)], 
(9) 
where r(Z, N) is the RMS charge radius of a nucleus with 
273 proton number Z and neutron number N. 


29 Figure 4 compares the experimental and calculated OES 


230 results for radii (A®) of the calcium (left panels) and cop- 
281 per (right panels) isotopes. The experimental data show that 
232 for the calcium isotopes [Figs. 4(a)-(e)] strong OES exists 
233 between N = 20 and 28 and that a reduction in the OES ap- 
284 pears for N > 28. Only RCHB theory could reproduce the 
trend of the experimental OES without KRR/EKRR correc- 
tions. However, the amplitude of the calculated OES was sig- 
nificantly less pronounced than that of the experimental data. 
Interestingly, after considering the KRR corrections, the cal- 
culated OES worsened for N < 28, particularly when the 
29 phase of the OES was opposite to that of the data. The A!/3- 
formula had no OES over the entire isotopic chain and the 
WS* model has a weak OES except at the N = 20 and 28 
shell closures. The OES in the HFB25 and HFB25* models 
were slightly higher. However, they were still weak compared 
with the data. Note that although OES can be obtained in 
the WS*, HFB25 and HFB25* models, the phases of the cal- 
culated OES are opposite to those of the experimental data. 
Considering the KRR method, the OES in these four mod- 
els increased, particularly for the WS* and HFB25* models 
for which the calculated OES were stronger than those of the 
data. However, the OES in these models were still opposite 
to those in the data. Therefore, although the KRR method 
improves the description of the charge radius to a large ex- 
tent, it was difficult to reproduce the observed OES. After 
considering the EKRR method, the experimental OES values 
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Fig. 3. (Color online) Comparison of the extrapolation ability of the KRR and EKRR methods for the neutron-rich region by considering six 
test sets with different extrapolation distances. The upper panels (a)-(e) show the RMS deviations of the KRR and EKRR methods. The lower 
panels (f)-(j) show the RMS deviations scaled to the corresponding RMS deviations for these five models without KRR/EKRR corrections. 
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Fig. 4. (Color online) Comparison of experimental and calculat- 
ed OES of the charge radii (A?) of the calcium (left panels) and 
copper (right panels) isotopes. The experimental data are shown 
as black squares. The calculation results of these five models are 
shown as grey solid circles, and the calculation results of the KR- 
R and EKRR models are shown as red triangles and blue crosses, 
respectively. 
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could be reproduced quite well, especially for the A!/* for- 
mula and RCHB theory for copper isotopes [Figs. 4(f)-(j)]. 
This situation is similar to that of the calcium isotopes. How- 
ever, the description of Cu isotopes is not as accurate as that 
of Ca isotopes when considering the EKRR corrections. The 
OES is overestimated in all these calculations for N < 33 
and N > 46. In addition, the phases of the OES between 
N = 38-40 were not well reproduced. However, the EKRR 
approach can improve the description of OES to a large ex- 
tent compared with the original theory. This indicates that 
after considering the odd-even effects, shell structures and 
many-body correlations, which are important for OES, can 
be learned well using an EKRR network. 


Similar to OES, abrupt kinks across the neutron shell clo- 
sures provide a particularly stringent test for nuclear theo- 
ry. In the present study, Pb and Sn isotopes were considered 
as examples for investigating the kinks across neutrons with 
N = 126 and 82 shell closures. Figure 5 compares the exper- 
imental and calculated differential mean-square charge radii. 
6(r2)N'N = (r2}N — (p2)N" for some, and even for Pb [Figs. 
5(a)-(e)] (relative to ?°8Pb, N’ = 126) and Sn [Figs. 5(f)-(j)] 
(relative to '8?Sn, N’ = 82). It can be observed that for Pb 
isotopes the RCHB theory can reproduce the kink at N = 126 
perfectly [Fig. 5(b)]. In the A!/® formula and HFB25 model, 
there is no kink [Figs. 5(a) and (c)]. The kink could be repro- 
duced using the WS* and HFB25* models, but with a slight 
overestimation [Figs. 5(d) and (e)]. The results obtained by 
considering the KRR and EKRR methods were similar. There 
are several interpretations of kinks [50, 82-85]. Our results 
indicate that kinks may not be connected to odd-even effects, 
such as pairing correlations. The well-reproduced kinks al- 
so provide a test of the proposed KRR/EKRR method. The 
kinks at N = 126 in all five models could be reproduced 
quite well, but the calculated differential mean-square charge 
radius at N = 132 was too large compared with the data. 
For the Sn isotopes, only the WS* and HFB25* models re- 
produced the kink at N = 82. However, the absolute values 
of the calculated 5(r?) from N = 74-78 are small compared 


377 
378 
379 
380 
381 
382 
383 
384 
385 
386 
387 
388 
389 
390 
391 
392 
393 
394 
395 


8(r’) 126/82, N [fm 


0.8 
0.4 
0.0 
-0.4 


124 128 
Neutron number 


132 74 76 78 80 82 84 


Neutron number 


Fig. 5. (Color online) Comparison of experimental and calculated d- 
ifferential mean-square charge radius Elroy NON = (r3 N — (r3 
for some even-even (a)-(e) Pb (relative to 2°°Pb, N’ = 126) and 
(f)-(j) Sn (relative to 18*Sn, N’ = 82) isotopes. The experimental 
data are shown as black squares. The results of these five model- 
s are shown as grey solid circles, and the calculation results of the 
KRR and EKRR models are shown as red triangles and blue crosses, 
respectively. 
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with the data, especially for the WS* model. After apply- 
ing the KRR/EKRR method, the results reproduced the data 
quite well. It also can be seen that the KRR/EKRR correction- 
s to the A!/3 formula and HFB25 model are inconspicuous. 
Therefore, the kink at N = 82 cannot be reproduced using the 
KRR/EKRR method. For the RCHB model, the differential 
mean-square charge radii calculated from N = 74 to-80 were 
improved, and a kink appeared, but was still slightly weaker 
compared with the data. 


V. SUMMARY 


In summary, the extended kernel ridge regression method 
with odd-even effects was adopted to improve the description 
of the nuclear charge radius by using five commonly used 
nuclear models. The hyperparameters of the KRR and EKRR 
methods for each model were determined using leave-one-out 
cross-validation. For each model, the resultant root-mean- 
square deviations of the 1014 nuclei with proton number 
Z > 8 could be significantly reduced to 0.009-0.013 fm after 
considering a modification with the EKRR method. The best 
among them was the RCHB model, with a root-mean-square 
deviation of 0.0092 fm, which is the best result for nuclear 
charge radii predictions using the machine learning approach 
as far as we know. The extrapolation abilities of the KRR and 
EKRR methods for the neutron-rich region were examined 
and it was found that after considering odd-even effects, the 
extrapolation power could be improved compared with that of 
the original KRR method. Strong odd-even staggering of nu- 
clear charge radii in Ca and Cu isotopes was investigated and 
reproduced quite well using the EKRR method. This indi- 
cates that after considering the odd-even effects, shell struc- 
tures and many-body correlations can be learned quite well 
using the EKRR network. Abrupt kinks across the neutron 
N = 126 and 82 shell closures were also investigated. 
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